Ultrahigh harmonics from laser-assisted ion-atom collisions 
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We present a theoretical analysis of high-order harmonic generation from ion-atom collisions 
in the presence of linearly polarized intense laser pulses. Photons with frequencies significantly 
higher than in standard atomic high-harmonic generation are emitted. These harmonics are due 
to two different mechanisms: (i) collisional electron capture and subsequent laser-driven transfer of 
an electron between projectile and target atom; (ii) reflection of a laser-driven electron from the 
projectile leading to recombination at the parent atom. 

PACS numbers: 42.65.Ky,34.50.Rk,32.80.Rm,34.70+c 



Over the last decades, a vast amount of work has been 
devoted to the study of ion-atom collisions |lj and atoms 
in intense laser fields However, the two areas were 

almost entirely separated. No experiments on ion-atom 
collisions in the presence of strong laser pulses have been 
carried out. The reported experiments on laser-assisted 
collisions involve one-photon processes and thermal 
collision energies. Also, theoretical descriptions |5| have 
mostly been limited to slow collisions and/or relatively 
weak fields. Recently, however, the theoretical works by 
M adsen et al. |6j and by Kirchner [7[ investigate fast col- 
lisions in the presence of a strong laser field. In Ref. 
excitation mechanisms are discussed, while Rcf. 7] fo- 
cuses on ionization and electron capture. In both cases, 
the presence of the field leads to a significant modifica- 
tion of the collision process. On the other hand, there has 
been no study on the question how typical strong-field 
processes in atoms, such as high-order harmonic gener- 
ation (HHG) y| and above-threshold ionization [jj, are 
modified due to the impact of an ion projectile. In HHG, 
a large number of incoming laser photons are converted 
into a single high-energy photon. HHG experiments are 
presently pursued with great effort [ToL Hi) since the pro- 
cess serves as a source of coherent XUV radiation and 
attosecond pulses. 

In the present work, we investigate HHG in laser- 
assisted ion-atom collision. We focus on impact veloc- 
ities such that the time-scales of nuclear and electronic 
motion are comparable, i.e., we have significant proba- 
bilities for collisional electron transfer from the target to 
the projectile. For sufficiently long laser pulse durations, 
the laser-driven electron effectively sees a large range of 
internuclear distances during one laser pulse. When the 
laser polarization axis is parallel to the ion impact veloc- 
ity, we show that this situation results in the generation 
of high harmonics with photon energies much higher than 
usually obtained in atomic HHG. 

The classical recollision model 01 describes atomic 
HHG as a sequence of strong-field ionization, acceler- 
ation of the electron in the laser field and recombina- 
tion with the core. Within this model, the maximum re- 



turn energy of the laser-driven electron is 3.17C/ p where 
Up = Eq I (4w 2 ) is the ponderomotive potential for a laser 
with field amplitude Eq and frequency to. The maximum 
energy of the emitted photons is then equal to 3.17U P +I p 
where J p is the atomic ionization potential. For ion-atom 
collisions, we show below how the interplay between colli- 
sional electron capture and laser-driven electron transfer 
between target and projectile leads to new mechanisms 
of HHG with cutoffs at significantly higher energies. 

We consider collisions of protons on hydrogen atoms for 
proton energies of 2 keV (impact velocity v — 0.283 a.u.). 
Due to the large impact momentum, the projectile trajec- 
tory is assumed to be classical and along a straight line. 
Furthermore, we use a two-dimensional model where all 
dynamics is restricted to the plane that contains the tar- 
get nucleus and the projectile. The internuclear vector 
is R(£) = (X, Z) = (b, vt) where b is the impact param- 
eter and v is the impact velocity. The interaction with 
the laser field E(f) = Eo(£) sm(ut + 5) is treated in the 
dipole approximation and in velocity gauge. The time- 
dependent Hamiltonian then reads (atomic units are used 
throughout) 



H(t) = ^- + p-A(t)+V(r t ) + V(r p ), 



(1) 



where A(t) 



f*~E(t')dt', r t = r + R(i)/2, and 



r p = r R(t)/2. For the electron-proton interaction 
V we use the softcore potential from Ref. For the 

laser field, we choose a wavelength of 800 nm and an 
intensity of lxlO 14 W/cm 2 since these parameters are 
readily available in experiment. We use a trapezoidal 
pulse shape with a three-cycle turn on and turn off. The 
total pulse duration is 16 optical cycles, which is equal 
to 42.7 fs. 

Initially, the electron is in the ground state of the target 
atom, i.e., localized around r t = 0. The time-dependent 
Schrodinger equation is then solved numerically by means 
of the split-operator method [1J| . The initial distance be- 
tween target and projectile is set to R = (b, —250 a.u.) 
so that the closest approach occurs at mid-pulse. 
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FIG. 1: Emission spectra for impact parameter b — 4 a.u. 
and laser phase 5 = 0. Full curve: polarization axis parallel to 
impact velocity. Dashed curve: polarization axis perpendicu- 
lar to impact velocity. The vertical lines indicate the position 
of cutoff energies for various mechanisms (see text). The inset 
shows the time-dependence of the laser field and the capture 
probability for the two polarization directions. 



The HHG spectrum is calculated according to 0] 



[a(t))e iuJt dt 



(2) 



where (a(t)) is the time-dependent dipole acceleration 
vector. 

Figure n shows the harmonic spectrum obtained for a 
collision with impact parameter 6 = 4 a.u. and phase 
S = of the laser. The two curves correspond to dif- 
ferent directions of the laser polarization. If the po- 
larization axis is perpendicular to the impact velocity, 
the emission spectrum has a form which is familiar from 
HHG in isolated atoms: a cutoff occurs at the photon en- 
ergy 3.17J7p + / p . Apparently, the passing projectile does 
not change the cutoff energy. If the polarization axis is 
parallel to the impact velocity, the result is strikingly 
different. We find an extension of HHG to frequencies 
reaching slightly beyond 32U P + I p . The region between 
3.17C/ p + I p and 8U P + I p appears like an extension of the 
atomic plateau with a steeper slope. (The significance 
of the values 8U P and 32U P will be explained below.) 
Furthermore, interesting hump structures appear in the 
spectrum around 3 a.u., 5.6 a.u., and 6.6 a.u. 

The inset of Fig. ^ shows the electric field E(t) and 
the capture probability P ca-P (t) which we define as the 
probability that the electron is found within a square of 
size 40 x 40 a.u. around the projectile. We see that 
the duration of the collisional capture process is much 
shorter than the pulse length since the ion-atom colli- 
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FIG. 2: Temporal analysis of harmonic emission for the same 
parameters as in Fig. The polarization axis is parallel to 
the impact velocity. Each spectrum describes the emission 
during one laser cycle and the value Z = vt refers to the 
middle of the respective time interval. The dotted horizontal 
lines indicate the respective level of log J = —10. The vertical 
lines indicate the same cutoff energies as in Fig. 



sion is essentially an attosecond process, cf. @ After 
capture at mid-pulse, P cap decreases due to ionization in 
the strong field. The capture probability is smaller for 
the case where the laser polarization is along the impact 
velocity, indicating that this geometry leads to more ion- 
ization during the ion-atom collision. 

Figure El gives a temporal analysis of the harmonic 
emission for the case that the laser is polarized paral- 
lel to the impact direction. Each spectrum in the fig- 
ure is obtained by Fourier transforming the dipole ac- 
celeration over one laser cycle. We find that the emis- 
sion below the atomic cutoff is nearly independent of 
time. The situation is not the same for the higher 
harmonics. Initially, there is no emission of ultrahigh 
harmonics (Z = —140.5 a.u.). Nevertheless, harmon- 
ics up to the highest frequency are produced already at 
Z = —109.3 a.u. and Z = —78.1 a.u., i.e. long before the 
actual ion-atom collision. Around Z = 0, the emission 
at the highest frequencies is weak. Instead, harmonics 
up to 8(7 P + I p are generated. At later times, emission 
at these energies drops and again, higher frequencies are 
produced. It is evident that the spectral structures at 
3 a.u., 5.6 a.u., and 6.6 a.u. in Fig.n mam ly arise from 
emission around Z = —46.8 a.u. and Z — —78.1 a.u. 
Since the generation of these harmonics is limited to al- 
most only one laser cycle, individual harmonics are not 
well resolved in those parts of the spectrum. The fact 
that ultrahigh harmonics are already generated at times 
when the approaching ion is farther than 100 a.u. from 
the target clearly shows that the presence of the ion in- 
fluences the dynamics of the laser-driven electrons when 
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FIG. 3: Panel (a): full curve, capture probability versus im- 
pact parameter for a laser field with phase 5 = 0; dashed 
curve, field- free capture probability. Panel (b): yield of ul- 
trahigh harmonics (full curve) and coherence parameter (dot- 
dashed curve) versus impact parameter. The polarization axis 
is parallel to the impact velocity. The insets show the same 
quantities as a function of the laser phase 5 for fixed impact 
parameter b — 4 a.u. 



they are far away from the target atom. 

As the next step in our analysis, we vary the impact 
parameter while keeping all other parameters constant. 
Figure OUa) displays the final capture probability P cap 
as a function of the impact parameter, with and with- 
out the presence of the laser. The capture probability is 
significantly reduced by strong-field ionization. The os- 
cillatory dependence on the impact parameter, however, 
is qualitatively unchanged. 

Since we are particularly interested in harmonic gener- 
ation well beyond the atomic cutoff, we define the yield 
of ultrahigh harmonics as the integrated quantity 



"Suh 



S(uj)duj. 



(3) 
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The solid line in Fig. Bib) displays the yield of ultrahigh 
harmonics as a function of the impact parameter. Its 
overall structure is a monotonic decrease. Thus, at first 
sight it exhibits no correlation with the capture prob- 
ability. A closer look, however, reveals that the "dips" 
around 6 — 2.5 a.u. and 6 = 4 a.u. coincide with extrema 
of the capture probability. For further investigation of 
this point we define the coherence parameter C accord- 
ing to C — PcupPt- Here, Pt is the final probability that 
the electron remains bound to the target atom, and is 
calculated analogous to P C ap- Small values of C indicate 
either ionization or that the electron is localized at one of 



the two nuclei. For large values of C, the electron state is 
a coherent superposition of target and projectile states. 
Figure Bib) shows that the maxima of C coincide with 
shoulders in S u h while the minima in C coincide with the 
dips in SxCn- 

We conclude that the generation of harmonics beyond 
the atomic cutoff originates from at least two different 
mechanisms. One mechanism is unrelated to the coher- 
ence parameter and gives rise to the monotonic decreas- 
ing background in Fig. Efb). A second mechanism de- 
pends on the strength of the coherence parameter and 
gives rise to the oscillations on top of the background. 

The insets in Fig. [3] show the dependence on the laser 
phase 5 for fixed impact parameter 6 = 4 a.u. As we 
vary the phase, we find only modest changes in the cap- 
ture probability, but a significant alteration of S u h. Note 
that phases between n and 2tt give a higher yield than 
phases between and ir. No clear correlation with the 
coherence parameter C is observed, indicating that the 
first mechanism dominates. This is consistent with the 
smallness of the coherence parameter at 6 = 4 a.u. 

As an explanation of our results, we propose the fol- 
lowing two mechanisms of HHG. In mechanism (i), one 
of two collision partners is ionized by the laser, the free 
electron is then accelerated in the time-dependent laser 
field, and finally the electron recombines with the other 
ion. Mechanism (ii) also begins with the creation and 
acceleration of a free electron. Instead of recombining 
at the other ion, however, the electron is elastically re- 
flected and is further accelerated in the laser field before 
it finally recombines with the core from which it was orig- 
inally ejected. 

Some years ago, mechanism (i) has been suggested 
as a new mechanism of HHG in stretched molecules 
[IH Iril Hif . and a cutoff at 8U P + I p was derived 
for the special internuclear distances R — (2n + V)tkol, 
n = 0, 1 . . ., where a = E /uj 2 is the classical quivering 
amplitude of the laser-driven electron. Until now, this 
cutoff has not been observed in experiment. One reason 
seems to be that rather large internuclear distances are 
required. Moreover, the effect occurs only for systems 
where the state of the electron prior to ionization can 
be described by a single-particle orbital that is coher- 
ently delocalized over both nuclei. For example, this is a 
valid description for the ground-state of an H^" molecular 
ion, but not for a neutral molecule at large internuclear 
distances. In an ion-atom collision, a coherent superpo- 
sition is realized if the coherence parameter C is appre- 
ciable. This explains the connection between C and S u h 
in Fig. Kb). 

Under the assumption of fixed nuclei, it is straightfor- 
ward to derive the maximum cutoff for both processes 
by inspection of the classical electronic equation of mo- 
tion, f = — Eosinwi. For an electron starting with zero 
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FIG. 4: Maximum kinetic electron energy versus internuclear 
distance for the recombination of a laser-driven electron in 
mechanisms (i) and (ii), see text (classical calculation for fixed 
nuclei) . 



velocity at r = at time to, we have 

f (t) = (E /w)(coswt — cos u>to). (4) 

Thus, the largest possible velocity equals 2Eq/u>, corre- 
sponding to an energy of 8U p . If an electron with this 
energy recombines at the other ion, it will generate a pho- 
ton with energy 8U P + I p . This explains the cutoff law 
for mechanism (i). If the electron is instead elastically 
backscattered at time t\, the velocity is thereafter given 
by 

r(t) = (Eo /lu) (cos ujt + cos u>to — 2coswii), t > t-y, (5) 

so that the maximum velocity is AEq /oj corresponding to 
an energy of 32t/ p . This maximum occurs for the case 
that to, t\, and t are times of zero electric field. It is 
straightforward to show that recollision with the maxi- 
mum energy is possible if the internuclear distance satis- 
fies the condition R = 3(2n + l)ira, n = 0, 1, . . . There- 
fore, at these distances we have a cutoff at 32U P + I p . 
For the present laser parameters, a = 16.5 a.u. and 
Up = 0.22 a.u. Although nonzero ion velocities may lead 
to a correction of the cutoff, we conclude that the gener- 
ation of ultrahigh harmonics (Figs. 1 1121) is well explained 
by the reflection mechanism. 

For mechanisms (i) and (ii), Fig. 01 shows the maxi- 
mum electron energy at recombination time as a function 
of the internuclear distance. These results are obtained 
by numerical solution of the classical equation of mo- 
tion for fixed nuclei. Both curves approach the atomic 
limit 3.17£/ p for R — > 0. The maximum energies are 
found at 32£7 p and 8U P , respectively, and this gives rise 
to the cutoffs mentioned above. Note also that already 
at R = ira, the reflection mechanism produces electron 
energies above 16U p . 

The phase dependence observed in Fig. [31 is consistent 
with the reflection scenario due to the asymmetry of the 



ion-atom system. Before the ion-atom collision, reflec- 
tion of electrons can occur only for electrons that are ac- 
celerated from the target towards the projectile. Under 
inversion of the field direction (phase shift by w) , these 
electrons are accelerated in the opposite direction where 
they cannot be backscattered. 

In summary, we have investigated ion-atom collisions 
in the presence of a strong laser field. High harmonics 
are generated with energies much larger than found in 
atomic HHG if the laser polarization axis is parallel to 
the direction of impact. We have proposed two distinct 
HHG mechanisms, both involving laser-driven transfer of 
electrons between the collision partners. The highest har- 
monics are produced by reflection of electrons from the 
projectile back to the target atom at times when target 
and projectile are far away from each other. A simplified 
classical description of the electron gives a cutoff energy 
at 32[/ p + I p for this process. In principle, the mecha- 
nisms of HHG that we have discussed should appear also 
in molecules. However, it seems difficult to combine the 
required large internuclear distances with precise control 
of the molecular orientation. In ion-atom collisions, all 
possible distances are sampled, and the orientation con- 
trol is straightforward. 
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